Early resource scarcity causes cortical astrocyte enlargement and sex-specific changes in the orbitofrontal cortex transcriptome in adult rats.

Astrocyte morphology affects function, including the regulation of glutamatergic signaling. This morphology changes dynamically in response to the environment. However, how early life manipulations alter adult cortical astrocyte morphology is underexplored. Our lab uses brief postnatal resource scarcity, the limited bedding and nesting (LBN) manipulation, in rats. We previously found that LBN augments maternal behaviors and promotes later resilience to adult addiction-related behaviors, reducing impulsivity, risky decision-making, and morphine self-administration. These behaviors rely on glutamatergic transmission in the medial orbitofrontal (mOFC) and medial prefrontal (mPFC) cortex. Here we tested whether LBN changed astrocyte morphology in the mOFC and mPFC of adult rats using a novel viral approach that, unlike traditional markers, fully labels astrocytes. Prior exposure to LBN causes an increase in the surface area and volume of astrocytes in the mOFC and mPFC of adult males and females relative to control-raised rats. We next used bulk RNA sequencing of OFC tissue to assess transcriptional changes that could increase astrocyte size in LBN rats. LBN caused mainly sex-specific changes in differentially expressed genes. Pathway analysis revealed that OFC glutamatergic signaling is altered by LBN in males and females, but the gene changes in that pathway differed across sex. This may represent a convergent sex difference where glutamatergic signaling, which affects astrocyte morphology, is altered by LBN via sex-specific mechanisms. Collectively, these studies highlight that astrocytes may be an important cell type that mediates the effect of early resource scarcity on adult brain function.

Early resource scarcity causes cortical astrocyte enlargement and sex-specific changes in the orbitofrontal cortex transcriptome in adult rats.

Astrocyte morphology affects function, including the regulation of glutamatergic signaling. This morphology changes dynamically in response to the environment. However, how early life manipulations alter adult cortical astrocyte morphology is underexplored. Our lab uses brief postnatal resource scarcity, the limited bedding and nesting (LBN) manipulation, in rats. We previously found that LBN promotes later resilience to adult addiction-related behaviors, reducing impulsivity, risky decision-making, and morphine self-administration. These behaviors rely on glutamatergic transmission in the medial orbitofrontal (mOFC) and medial prefrontal (mPFC) cortex. Here we tested whether LBN changed astrocyte morphology in the mOFC and mPFC of adult rats using a novel viral approach that, unlike traditional markers, fully labels astrocytes. Prior exposure to LBN causes an increase in the surface area and volume of astrocytes in the mOFC and mPFC of adult males and females relative to control-raised rats. We next used bulk RNA sequencing of OFC tissue to assess transcriptional changes that could increase astrocyte size in LBN rats. LBN caused mainly sex-specific changes in differentially expressed genes. However, Park7, which encodes for the protein DJ-1 that alters astrocyte morphology, was increased by LBN across sex. Pathway analysis revealed that OFC glutamatergic signaling is altered by LBN in males and females, but the gene changes in that pathway differed across sex. This may represent a convergent sex difference where glutamatergic signaling, which affects astrocyte morphology, is altered by LBN via sex-specific mechanisms. Collectively, these studies highlight that astrocytes may be an important cell type that mediates the effect of early resource scarcity on adult brain function.

Paternal morphine exposure enhances morphine self-administration and induces region-specific neural adaptations in reward-related brain regions of male offspring.

Background: A growing body of preclinical studies report that preconceptional experiences can have a profound and long-lasting impact on adult offspring behavior and physiology. However, less is known about paternal drug exposure and its effects on reward sensitivity in the next generation. Methods: Adult male rats self-administered morphine for 65 days; controls received saline. Sires were bred to drug-naïve dams to produce first-generation (F1) offspring. Morphine, cocaine, and nicotine self-administration were measured in adult F1 progeny. Molecular correlates of addiction-like behaviors were measured in reward-related brain regions of drug naïve F1 offspring. Results: Male, but not female offspring produced by morphine-exposed sires exhibited dose-dependent increased morphine self-administration and increased motivation to earn morphine infusions under a progressive ratio schedule of reinforcement. This phenotype was drug-specific as self-administration of cocaine, nicotine, and sucrose were not altered by paternal morphine history. The male offspring of morphine-exposed sires also had increased expression of mu-opioid receptors in the ventral tegmental area but not in the nucleus accumbens. Conclusions: Paternal morphine exposure increased morphine addiction-like behavioral vulnerability in male but not female progeny. This phenotype is likely driven by long-lasting neural adaptations within the reward neural brain pathways.

A systems biology analysis of protein-protein interaction of digestive disorders and Covid-19 virus based on comprehensive gene information.

Aim: Analysis of networks of digestive disorder and their relationship with Covid-19 based on systems biology methods, evaluation similarity, and usefulness of networks to give a new treatment approach. Background: Digestive disorders are typically complex diseases associated with high treatment costs. They are related to the immune system and inflammation. With the outbreak of Covid-19, this disease was shown to have signs like diarrhea. Some signs of Covid-19 are similar to those of digestive disorders, like IBD and diarrhea. Both of them are accompanied by inflammation and induce disorders in the digestive system. Methods: DisGeNET and STRING databases were sources of disease genes and constructing networks and were used to construct the network of digestive diseases and Covid-19. Three plugins of Cytoscape software, namely ClusterONE, ClueGO, and CluePedia, were used to analyze cluster networks and enrichment pathways. To describe the interaction of proteins, information from KEGG pathway and Reactome was used. Results: According to the results, IBD, gastritis, and diarrhea have common pathways. The CXCL8, IL-6, IL-1ß, TNF-α, TLR4, and MBL2 molecules were identified as inflammatory molecules in all networks. Conclusion: It seems that detecting genes and pathways can be useful in applying new approaches for treating these diseases.

Transcriptomics in the nucleus accumbens shell reveal sex- and reinforcer-specific signatures associated with morphine and sucrose craving.

Incubation of craving is a well-documented phenomenon referring to the intensification of drug craving over extended abstinence. The neural adaptations that occur during forced abstinence following chronic drug taking have been a topic of intense study. However, little is known about the transcriptomic changes occurring throughout this window of time. To define gene expression changes associated with morphine consumption and extended abstinence, male and female rats underwent 10 days of morphine self-administration. Separate drug-naive rats self-administered sucrose in order to compare opioid-induced changes from those associated with natural, non-drug rewards. After one or 30 days of forced abstinence, rats were tested for craving, or nucleus accumbens shell tissue was dissected for RNA sequencing. Morphine consumption was predictive of drug seeking after extended (30 days) but not brief (1 day) abstinence in both sexes. Extended abstinence was also associated with robust sex- and reinforcer-specific changes in gene expression, suggesting sex differences underlying incubation of morphine and sucrose seeking respectively. Importantly, these changes in gene expression occurred without re-exposure to drug-paired cues, indicating that chronic morphine causes long-lasting changes in gene expression that prime the system for increased craving. These findings lay the groundwork for identifying specific therapeutic targets for curbing opioid craving without impacting the natural reward system in males and females.

Effects of early life adversity on male reproductive behavior and the medial preoptic area transcriptome.

Early life adversity can alter reproductive development in humans, changing the timing of pubertal onset and sexual activity. One common form of early adversity is limited access to resources. This adversity can be modeled in rats using the limited bedding/nesting model (LBN), in which dams and pups are placed in a low resource environment from pups' postnatal days 2-9. Our laboratory previously found that adult male rats raised in LBN conditions have elevated levels of plasma estradiol compared to control males. In females, LBN had no effect on plasma hormone levels, pubertal timing, or estrous cycle duration. Estradiol mediates male reproductive behaviors. Thus, here we compared reproductive behaviors in adult males exposed to LBN vs. control housing. LBN males acquired the suite of reproductive behaviors (mounts, intromissions, and ejaculations) more quickly than their control counterparts over 3 weeks of testing. However, there was no effect of LBN in males on puberty onset or masculinization of certain brain regions, suggesting LBN effects on estradiol and reproductive behaviors manifest after puberty. In male and female rats, we next used RNA sequencing to characterize LBN-induced transcriptional changes in the medial preoptic area (mPOA), which underlies male reproductive behaviors. LBN produced sex-specific alterations in gene expression, with many transcripts showing changes in opposite directions. Numerous transcripts altered by LBN in males are regulated by estradiol, linking hormonal changes to molecular changes in the mPOA. Pathway analysis revealed that LBN induced changes in neurosignaling and immune signaling in males and females, respectively. Collectively, these studies reveal novel neurobiological mechanisms by which early life adversity can alter reproductive strategies.

Chronic paternal morphine exposure increases sensitivity to morphine-derived pain relief in male progeny.

Parental history of opioid exposure is seldom considered when prescribing opioids for pain relief. To explore whether parental opioid exposure may affect sensitivity to morphine in offspring, we developed a "rat pain scale" with high-speed imaging, machine learning, and mathematical modeling in a multigenerational model of paternal morphine self-administration. We find that the most commonly used tool to measure mechanical sensitivity in rodents, the von Frey hair, is not painful in rats during baseline conditions. We also find that male progeny of morphine-treated sires had no baseline changes in mechanical pain sensitivity but were more sensitive to the pain-relieving effects of morphine. Using RNA sequencing across pain-relevant brain regions, we identify gene expression changes within the regulator of G protein signaling family of proteins that may underlie this multigenerational phenotype. Together, this rat pain scale revealed that paternal opioid exposure increases sensitivity to morphine's pain-relieving effects in male offspring.

Novel insight into pancreatic adenocarcinoma pathogenesis using liquid association analysis.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy associated with a poor prognosis. High-throughput disease-related-gene expression data provide valuable information on gene interaction, which consequently lead to deeper insight about pathogenesis. The co-expression analysis is a common approach that is used to investigate gene interaction. However, such an approach solely is inadequate to reveal the complexity of the gene interaction. The three-way interaction model is known as a novel approach applied to decode the complex relationship between genes. METHODS: In the current study, the liquid association method was used to capture the statistically significant triplets involved in the PDAC pathogenesis. Subsequently, gene set enrichment and gene regulatory network analyses were performed to trace the biological relevance of the statistically significant triplets. RESULTS: The results of the current study suggest that "response to estradiol" and "Regulation of T-cell proliferation" are two critical biological processes that may be associated with the PDAC pathogenesis. Additionally, we introduced six switch genes, namely Lamc2, Klk1, Nqo1, Aox1, Tspan1, and Cxcl12, which might be involved in PDAC triggering. CONCLUSION: In the current study, for the first time, the critical genes and pathways involved in the PDAC pathogenesis were investigated using the three-way interaction approach. As a result, two critical biological processes, as well as six potential biomarkers, were suggested that might be involved in the PDAC triggering. Surprisingly, strong evidence for the biological relevance of our results can be found in the literature.

IrGO: Iranian traditional medicine General Ontology and knowledge base.

BACKGROUND: Iranian traditional medicine, also known as Persian Medicine, is a holistic school of medicine with a long prolific history. It describes numerous concepts and the relationships between them. However, no unified language system has been proposed for the concepts of this medicine up to the present time. Considering the extensive terminology in the numerous textbooks written by the scholars over centuries, comprehending the totality of concepts is obviously a very challenging task. To resolve this issue, overcome the obstacles, and code the concepts in a reusable manner, constructing an ontology of the concepts of Iranian traditional medicine seems a necessity. CONSTRUCTION AND CONTENT: Makhzan al-Advieh, an encyclopedia of materia medica compiled by Mohammad Hossein Aghili Khorasani, was selected as the resource to create an ontology of the concepts used to describe medicinal substances. The steps followed to accomplish this task included (1) compiling the list of classes via examination of textbooks, and text mining the resource followed by manual review to ensure comprehensiveness of extracted terms; (2) arranging the classes in a taxonomy; (3) determining object and data properties; (4) specifying annotation properties including ID, labels (English and Persian), alternative terms, and definitions (English and Persian); (5) ontology evaluation. The ontology was created using Protégé with adherence to the principles of ontology development provided by the Open Biological and Biomedical Ontology (OBO) foundry. UTILITY AND DISCUSSION: The ontology was finalized with inclusion of 3521 classes, 15 properties, and 20,903 axioms in the Iranian traditional medicine General Ontology (IrGO) database, freely available at http://ir-go.net/ . An indented list and an interactive graph view using WebVOWL were used to visualize the ontology. All classes were linked to their instances in UNaProd database to create a knowledge base of ITM materia medica. CONCLUSION: We constructed an ontology-based knowledge base of ITM concepts in the domain of materia medica to help offer a shared and common understanding of this concept, enable reuse of the knowledge, and make the assumptions explicit. This ontology will aid Persian medicine practitioners in clinical decision-making to select drugs. Extending IrGO will bridge the gap between traditional and conventional schools of medicine, helping guide future research in the process of drug discovery.

Identification of Celecoxib-Targeted Proteins Using Label-Free Thermal Proteome Profiling on Rat Hippocampus.

Celecoxib, or Celebrex, a nonsteroidal anti-inflammatory drug, is one of the most common medicines for treating inflammatory diseases. Recently, it has been shown that celecoxib is associated with implications in complex diseases, such as Alzheimer disease and cancer as well as with cardiovascular risk assessment and toxicity, suggesting that celecoxib may affect multiple unknown targets. In this project, we detected targets of celecoxib within the nervous system using a label-free thermal proteome profiling method. First, proteins of the rat hippocampus were treated with multiple drug concentrations and temperatures. Next, we separated the soluble proteins from the denatured and sedimented total protein load by ultracentrifugation. Subsequently, the soluble proteins were analyzed by nano-liquid chromatography tandem mass spectrometry to determine the identity of the celecoxib-targeted proteins based on structural changes by thermal stability variation of targeted proteins toward higher solubility in the higher temperatures. In the analysis of the soluble protein extract at 67°C, 44 proteins were uniquely detected in drug-treated samples out of all 478 identified proteins at this temperature. Ras-associated binding protein 4a, 1 out of these 44 proteins, has previously been reported as one of the celecoxib off targets in the rat central nervous system. Furthermore, we provide more molecular details through biomedical enrichment analysis to explore the potential role of all detected proteins in the biologic systems. We show that the determined proteins play a role in the signaling pathways related to neurodegenerative disease-and cancer pathways. Finally, we fill out molecular supporting evidence for using celecoxib toward the drug-repurposing approach by exploring drug targets. SIGNIFICANCE STATEMENT: This study determined 44 off-target proteins of celecoxib, a nonsteroidal anti-inflammatory and one of the most common medicines for treating inflammatory diseases. It shows that these proteins play a role in the signaling pathways related to neurodegenerative disease and cancer pathways. Finally, the study provides molecular supporting evidence for using celecoxib toward the drug-repurposing approach by exploring drug targets.

Early life adversity promotes resilience to opioid addiction-related phenotypes in male rats and sex-specific transcriptional changes.

Experiencing some early life adversity can have an "inoculating" effect that promotes resilience in adulthood. However, the mechanisms underlying stress inoculation are unknown, and animal models are lacking. Here we used the limited bedding and nesting (LBN) model of adversity to evaluate stress inoculation of addiction-related phenotypes. In LBN, pups from postnatal days 2 to 9 and their dams were exposed to a low-resource environment. In adulthood, they were tested for addiction-like phenotypes and compared to rats raised in standard housing conditions. High levels of impulsivity are associated with substance abuse, but in males, LBN reduced impulsive choice compared to controls. LBN males also self-administered less morphine and had a lower breakpoint on a progressive ratio reinforcement schedule than controls. These effects of LBN on addiction-related behaviors were not found in females. Because the nucleus accumbens (NAc) mediates these behaviors, we tested whether LBN altered NAc physiology in drug-naïve and morphine-exposed rats. LBN reduced the frequency of spontaneous excitatory postsynaptic currents in males, but a similar effect was not observed in females. Only in males did LBN prevent a morphine-induced increase in the AMPA/NMDA ratio. RNA sequencing was performed to delineate the molecular signature in the NAc associated with LBN-derived phenotypes. LBN produced sex-specific changes in transcription, including in genes related to glutamate transmission. Collectively, these studies reveal that LBN causes a male-specific stress inoculation effect against addiction-related phenotypes. Identifying factors that promote resilience to addiction may reveal novel treatment options for patients.

Can We Assume the Gene Expression Profile as a Proxy for Signaling Network Activity?

Studying relationships among gene products by expression profile analysis is a common approach in systems biology. Many studies have generalized the outcomes to the different levels of central dogma information flow and assumed a correlation of transcript and protein expression levels. However, the relation between the various types of interaction (i.e., activation and inhibition) of gene products to their expression profiles has not been widely studied. In fact, looking for any perturbation according to differentially expressed genes is the common approach, while analyzing the effects of altered expression on the activity of signaling pathways is often ignored. In this study, we examine whether significant changes in gene expression necessarily lead to dysregulated signaling pathways. Using four commonly used and comprehensive databases, we extracted all relevant gene expression data and all relationships among directly linked gene pairs. We aimed to evaluate the ratio of coherency or sign consistency between the expression level as well as the causal relationships among the gene pairs. Through a comparison with random unconnected gene pairs, we illustrate that the signaling network is incoherent, and inconsistent with the recorded expression profile. Finally, we demonstrate that, to infer perturbed signaling pathways, we need to consider the type of relationships in addition to gene-product expression data, especially at the transcript level. We assert that identifying enriched biological processes via differentially expressed genes is limited when attempting to infer dysregulated pathways.

UNaProd: A Universal Natural Product Database for Materia Medica of Iranian Traditional Medicine.

BACKGROUND: Iranian traditional medicine (ITM) is a holistic medical system that uses a wide range of medicinal substances to treat disease. Reorganization and standardization of the data on ITM concepts is a necessity for optimal use of this rich source. In an initial step towards this goal, we created a database of ITM materia medica. Main Body. Primarily based on Makhzan al-Advieh, which is the most recent encyclopedia of materia medica in ITM with the largest number of monographs, a database of natural medicinal substances was created using both text mining methods and manual editing. UNaProd, a Universal Natural Product database for materia medica of ITM, is currently host to 2696 monographs, from herbal to animal to mineral compounds in 16 diverse attributes such as origin and scientific name. Currently, systems biology, and more precisely systems medicine and pharmacology, can be an aid in providing rationalizations for many traditional medicines and elucidating a great deal of knowledge they can offer to guide future research in medicine. CONCLUSIONS: A database of materia medica is a stepping stone in creating a systems pharmacology platform of ITM that encompasses the relationships between the drugs, their targets, and diseases. UNaProd is hyperlinked to IrGO and CMAUP databases for Mizaj and molecular features, respectively, and it is freely available at http://jafarilab.com/unaprod/.

Showing NAFLD, as a key connector disease between Alzheimer's disease and diabetes via analysis of systems biology.

AIM: This study was designed to perform network analysis of Alzheimers҆ disease and diabetes and to find their correlation with each other and other diseases/pathways. BACKGROUND: Alzheimer's disease (AD) as a neurodegenerative disease and diabetes as a metabolic disease are two major health problems in the recent years. The recent studies have reported their correlation and same spreading pathways of these two diseases together, but details of this relation are not well known yet at molecular level.. METHODS: In thermal proteome profiling (TPP) technique, after treatment of the extracted proteins by heat and drug concentration, the resulting proteins were analyzed by mass spectrometry. Enrichment analysis of these proteins led to development of AD and diabetes. First, corresponding genes for each disease were extracted from DisGeNET database and then, protein-protein interaction network was constructed for each of them using the search tool for retrieval of interacting genes and proteins (STRING). After analyzing these networks, hub-bottleneck nodes of networks were evaluated. Also, common nodes between two networks were extracted and used for further analysis. RESULTS: High correlation was found between AD and diabetes based on the existence of 40 common genes. Results of analyses revealed 14 genes in AD and 12 genes in diabetes as hub-bottleneck 7 of which were common including caspase 3 (CASP3), insulin-like growth factor 1 (IGF1), catalase (CAT), tumor necrosis factor (TNF), leptin (LEP), vascular endothelial growth factor A (VEGFA), and interleukin 6 ( IL-6). CONCLUSION: Our results revealed a direct correlation between AD and diabetes and also a correlation between these two diseases and non-alcoholic fatty liver disease (NAFLD), suggesting that a small change in each of these three diseases can lead to development of any other diseases in the patients. Also, the enrichments exhibited the existence of common pathways between AD, diabetes, NAFLD, and male infertility.

The developmental origins of sex-biased expression in cardiac development.

BACKGROUND: Expression patterns between males and females vary in every adult tissue, even in organs with no conspicuous dimorphisms such as the heart. While studies of male and female differences have traditionally focused on the influence of sex hormones, these do not account for all the differences at the molecular and epigenetic levels. We previously reported that a substantial number of genes were differentially expressed in male and female mouse embryonic stem (ES) cells and revealed dose-dependent enhancer activity in response to Prdm14, a key pluripotency factor expressed more highly in female ES cells. In this work, we investigated the role of Prdm14 in establishing sex-specific gene expression networks. We surveyed the sex-specific landscape in early embryogenesis with special reference to cardiac development. We generated sex-specific co-expression networks from mouse ES cells, examined the presence of sex-specific chromatin domains, and analyzed previously published datasets from different developmental time points to characterize how sex-biased gene expression waxes and wanes to evaluate whether sex-biased networks are detectable throughout heart development. RESULTS: We performed ChIP-seq on male and female mouse ES cells to determine differences in chromatin status. Our study reveals sex-biased histone modifications, underscoring the potential for the sex chromosome complement to prime the genome differently in early development with consequences for later expression biases. Upon differentiation of ES cells to cardiac precursors, we found sex-biased expression of key transcription and epigenetic factors, some of which persisted from the undifferentiated state. Using network analyses, we also found that Prdm14 plays a prominent role in regulating a subset of dimorphic expression patterns. To determine whether sex-biased expression is present throughout cardiogenesis, we re-analyzed data from two published studies that sampled the transcriptomes of mouse hearts from 8.5 days post-coitum embryos to neonates and adults. We found sex-biased expression at every stage in heart development, and interestingly, identified a subset of genes that exhibit the same bias across multiple cardiogenic stages. CONCLUSIONS: Overall, our results support the existence of sexually dimorphic gene expression profiles and regulatory networks at every stage of cardiac development, some of which may be established in early embryogenesis and epigenetically perpetuated.

Protein-protein interaction analysis of Alzheimer`s disease and NAFLD based on systems biology methods unhide common ancestor pathways.

AIM: Analysis reconstruction networks from two diseases, NAFLD and Alzheimer`s diseases and their relationship based on systems biology methods. BACKGROUND: NAFLD and Alzheimer`s diseases are two complex diseases, with progressive prevalence and high cost for countries. There are some reports on relation and same spreading pathways of these two diseases. In addition, they have some similar risk factors, exclusively lifestyle such as feeding, exercises and so on. Therefore, systems biology approach can help to discover their relationship. METHODS: DisGeNET and STRING databases were sources of disease genes and constructing networks. Three plugins of Cytoscape software, including ClusterONE, ClueGO and CluePedia, were used to analyze and cluster networks and enrichment of pathways. An R package used to define best centrality method. Finally, based on degree and Betweenness, hubs and bottleneck nodes were defined. RESULTS: Common genes between NAFLD and Alzheimer`s disease were 190 genes that used construct a network with STRING database. The resulting network contained 182 nodes and 2591 edges and comprises from four clusters. Enrichment of these clusters separately lead to carbohydrate metabolism, long chain fatty acid and regulation of JAK-STAT and IL-17 signaling pathways, respectively. Also seven genes selected as hub-bottleneck include: IL6, AKT1, TP53, TNF, JUN, VEGFA and PPARG. Enrichment of these proteins and their first neighbors in network by OMIM database lead to diabetes and obesity as ancestors of NAFLD and AD. CONCLUSION: Systems biology methods, specifically PPI networks, can be useful for analyzing complicated related diseases. Finding Hub and bottleneck proteins should be the goal of drug designing and introducing disease markers.

A systems biology analysis protein-protein interaction of NASH and IBD based on comprehensive gene information.

AIM: Analysis reconstruction networks from two diseases, IBD and NASH and their relationship, based on systems biology methods. BACKGROUND: IBD and NASH are two complex diseases, with progressive prevalence and high cost for countries. There are some reports on co-existence of these two diseases. In addition, they have some similar risk factors such as age, obesity, and insulin resistance. Therefore, systems biology approach can help to discover their relationship. METHODS: DisGeNET and STRING databases were sources of disease genes and constructing networks. Three plugins of Cytoscape software, including ClusterONE, ClueGO and CluePedia, were used to analyze and cluster networks and enrichment of pathways. Based on degree and Betweenness, hubs and bottleneck nodes were defined. RESULTS: Common genes between IBD and NASH construct a network with 99 nodes. Common genes between IBD and NASH were extracted and imported to STRING database to construct PPI network. The resulting network contained 99 nodes and 333 edges. Five genes were selected as hubs: JAK2, TLR2, TP53, TLR4 and STAT3 and five genes were selected as bottleneck including: JAK2, TP53, AGT, CYP3A4 and TLR4. These genes were hubs in analysis network that was constructed from hubs of NASH and IBD networks. CONCLUSION: Systems biology methods, specifically PPI networks, can be useful for analyzing complicated related diseases. Finding Hub and bottleneck proteins should be the goal of drug designing and introducing disease markers.